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Abstract 

Bessel beams are studied within the general framework of quantum optics. The two modes of the 
electromagnetic field are quantized and the basic dynamical operators are identified. The algebra 
of these operators is analyzed in detail; it is shown that the operators that are usually associated to 
linear momentum, orbital angular momentum and spin do not satisfy the algebra of the translation 
and rotation group, and that this algebra is not closed. Some physical consequences of these results 
are examined. 

PACS numbers: 42.50.Vk, 32.80.Lg 
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I. INTRODUCTION 



The increasing use of light to control the motion of atomic systems and microparticles 
has renewed the interest in the mechanical properties of electromagnetic (EM) beams. The 
realization that light carries energy as well as linear and angular momenta was essential 
in the development of the classical theory of electromagnetism; nowadays, the interchange 
of such quantities with matter is a well established fact. In most cases, the interaction of 
light with matter can be satisfactorily described by taking the EM field as a superposition 
of idealized plane waves; in this simple picture, each normal mode carries momenta in the 
direction of the propagation vector, and the angular momentum is directly related to the 
states of polarization. 

However, the decomposition of angular momentum into an orbital and spin part has some 
ambiguities within the framework of quantum optics, as was recognized in classical papers 
by Darwin [lj], Humblet Q], de Broglie j^l, and many others. Recent interest in defining 
angular momentum in optics can be traced back to the works of Lenstra and Mandel |4|, 
Allen et al jjj, and van Enk and Nienhuis jfj, Q|. Lenstra and Mandel considered periodic 
boundary conditions that limit the isotropy of the EM field and thus affect the angular 
momentum. Allen et al showed that angular momentum can indeed be decomposed into 
orbital and spin parts for Laguerre-Gaussian modes, i. e., paraxial elementary waves with 
cylindrical symmetry. Finally, van Enk and Nienhuis [t| studied this decomposition and 
showed that neither the spin nor the orbital quantum operators of the electromagnetic field 
satisfy the commutation relations of angular momentum; they made the assumption that 
the total angular momentum obtained from Noether theorem is given directly as a sum of 
these operators; however, this is not the case in many situations where boundary conditions 
give rise to important surface effects. 

In general, orbital angular momentum (OAM) is identified with the part of the total 
angular momentum that depends explicitly on position (with respect to an origin of coor- 
dinates). When OAM is evaluated with respect to the axis of symmetry, it turns out that 
Laguerre-Gaussian modes do carry OAM in integer multiples of h jjj. However, the formal 
decomposition into OAM and spin angular momentum does not appear to be natural beyond 
the paraxial approximation Nevertheless, recent experiments have shown that angular 
momentum do indeed possess orbital and spin parts Q, H, 12, the former having 
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an intrinsic and extrinsic nature with direct physical consequences |15|- 

Bessel beams have interesting properties that make them especially attractive: they prop- 
agate with an intensity pattern invariant along a given axis and carry angular momentum 
along that axis. Experimental realizations of such beams and their mechanical effects on 



atoms and microparticles are the subject of many current investigations kJ. The purpose 
of the present paper is to study the mechanical properties of Bessel beams within the general 
formalism of quantum optics. In particular, we study the standard decomposition of the 
total angular momentum of these beams into orbital and spin parts, and obtain explicit 
expressions for these observables as quantum operators. It turns out, however, that they 
do not satisfy the usual commutation relations for angular momentum vectors and that, 
contrary to van Enk and Nienhuis (j| results, the algebra does not even close; we argue that 
this discrepancy is due to the boundary conditions. We also show that the superposition of 
Bessel modes, as well as their polarization states, can be characterized by a set of opera- 
tors that appears naturally within our formulation and can be measured in principle. Some 
predictions that could eventually be tested experimentally are discussed in the conclusions. 

The organization of the article is as follows. In section II the electromagnetic Bessel 
beams are expressed in terms of Hertz potentials. In Section III the fields are quantized 
and explicit expressions for the operators are given, relating them to their main mechanical 
properties. Section IV is devoted to the study of the algebraic properties of the operators that 
are usually identified with orbital and spin angular momentum. Finally, a brief discussion 
of the results and their experimental implications is presented in Section V. Some useful 
formulas are given in Appendix A, and expressions relating Bessel modes with plane and 
spherical EM modes are given in Appendix B. 

II. ELECTROMAGNETIC BESSEL MODES 

An electromagnetic field with cylindrical symmetry can be conveniently described in the 
terms of Hertz potentials X and G2 }rfl ]. In cylindrical coordinates {p, <p, z}, the electro- 
magnetic potentials are given by 

• = -!*, (i) 
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. p dcj) ' dp ' cdt 
and satisfy the Lorentz gauge condition. Then, the electric field E is: 
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E <t> = -7T^ i + ^r 2, (4) 



p dzdcj) cdtdp 
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and the magnetic field B is: 



* = + 936., (5) 



1 a 2 a 2 

B p = -^- i + 7Tir 2, (6) 



p cdtdcj) dzdp 
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cdtdp^ 1 ^ pdzdc^ 2 '' ^ 



a 2 a 2 

^ = -3^© 2 + 7^©2. (8) 



c 2 <9t 2 <9z 2 

Both Hertz potentials Gj, (i = 1,2), satisfy the equations: 

d 2 „ 1 d ( d „ \ 1 d 2 „ d 2 



c 2 <9t 2 pdp\ dp ) p 2 d(f) 2 dz 2 
Any solution of this equation that is regular at the origin can be written as a linear 
combination of the functions 

0i = CiJ m (k±p) exp{— iut + ik z z + imcj)}, (10) 



where J m is the Bessel function of order m, Cj are constants, and A;^ = a/ (y/c) 2 — A; 2 . Bessel 
functions form a complete orthogonal basis as follows from Eq. (|73|) in Appendix A. 

An electromagnetic mode is associated to each Hertz potential, Gi and G2, via Eqs. (0|H|), 
giving rise to transverse magnetic and electric modes respectively. In the following we will 
occasionally denote them by the superscripts TM and TE. 

It is convenient to make a further gauge transformation to the transverse gauge with 
$ = 0. This can be achieved with the transformations $ — > $ — dA/cdt and A — > A + VA, 
and taking A = (k z c/u)<d 1 . 



The electromagnetic vectors and potential can be decomposed in terms of their basic 
modes. First, we define the vectors: 



M(r, t;K) = ^- ■^J m {k±p)e p + % J' m (k ± p)e rt> 



-iujt+imtfr+ikz z 



and 



m 



N(r,t;K) = U' m (k ± p)e p - ■^J m (k ± p)e 4 



k± 

+ -rrJ m (k±p)e z 

K 7 , 



-iwt+imrfi+ikz z 



(11) 



(12) 



Here and in the following, the set of quantum numbers {k±, m, k z } are denoted by the generic 
symbol K whenever no confusion can arise. 

Accordingly, we obtain the following forms for the modes of the electromagnetic potential: 



A(™\r,t;K) = -£j? M \k ± , k z )N(v, t; K) 



iuj 



(13) 



and 



A< TS >(r, t; K) = --£^ E \k ± , k z )M(v, t; K), (14) 

where the functions £m(k±,k z ) (i = TM,TE) refer to the amplitudes of the transverse 
electric and magnetic modes. 

The electric field for each mode is now given by: 

E«(r,t; K) = — A w (r,t; K), 



and the magnetic field is: 

B(™\v,t;K) = ££ M \k ± ,k z )M(r,t;K), 
B^ E \r,t;K) = £% E \k ± , fc,)N(r, t; K). 

Finally, we notice that the vector N can also be written in the form: 

N{v,t-K) = - l -\j m+1 (k ± p)e^ m+1 ^ - J m -i(k ± p)e^ m ^e 



-iu)t+ik z z 



(15) 

(16) 
(17) 



with 



+-±Jm(k±p)e tm *e 3 , 

K z 



e± = ei ± ie 2 , 



(18) 
(19) 
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and that ck z ~M = lu~N x e 3 . Some useful formulas involving the vectors M and N are given 
in the Appendices. 

Notice also that the electromagnetic field described by the above expressions is purely- 
transverse, in the sense that V • E = 0. 

For the sake of comparison, we have included an appendix with the expansion of the 
Bessel modes in terms of the more common plane and spherical waves. 



III. QUANTIZATION AND DYNAMICAL VARIABLES. 

In the formalism of the previous section, the electromagnetic field is described in terms of 
two independent sets of modes. This representation has the advantage that the field can be 
quantized without further complications, even though it is a vector field with an additional 
freedom of gauge. The field operator A(r, t) takes the explicit form: 

00 poo poo 

A M) = E E / dk ± / dk z [a%(k z ,k ± )A^( r ,t;K) 

i=l,2m=-oo^° 

+ a^(k z ,k ± )A^*(v,t;K)}, (20) 
where the annihilation and creation operators satisfy the usual commutation relations 

[aS(k ± , k z ), a%)\k' ± , k' z )] = 5^5 m , m ,5(k ± - k' ± )5(k z - k' z ), (21) 

the index i referring to the two modes of the electromagnetic field, that is, the TM(i = 1) 
and TEii = 2) modes. 

The modes should be so normalized that each photon of frequency 

uj = c\J k 2 z + k\ 

carries an energy ftw given by 

^(^) = ^ / [|E(r,t;i^)| 2 + |B(r,t;^T)| 2 ]^ . (22) 

This condition can be satisfied provided the integral is well defined. For Bessel modes, with 
the volume of integration taken as the whole space, the normalization condition must be 
generalized to: 

J [E^*(r,t;K) -E w (r,t; K') + B w *(r,t; K) ■ B (i) (r, t; K')] dV 



= hw5 m ^5(k± - k' ± )5(k z - k' z ) , (23) 

which is equivalent to choosing amplitudes £m E \k±,k z ) = £m (k±,k z ) = k z c\Jhk±/2ixLLi 
for each mode. 

Defining now the generalized number operator. 

= 2 {^m^m + a m^m^5 (24) 

the quantum energy operator takes the form: 

£ = hJ2 [ dk ± dk z u N®(k±,k g ), (25) 

as it should be. 

For time independent states, the expectation value of the energy operator, integrated 
over a certain volume, will remain constant as long as the total energy flux over a surface 
enclosing that volume is zero. If the volume is the whole space, this condition is satisfied 
provided the field is localized, that is, its expectation value decays to zero sufficiently fast 
at infinity. 

We now turn our attention to other dynamical variables. The general expression for the 
momentum operator is 

P(t) = ^- / [E(r,t) xB(r,t)-B(r,t) x E(r, t)]dV. (26) 
Stic J L J 

For the Bessel modes under consideration, it takes the form: 

P = hJ2 I dk ± dk z [ik ± a { Sid%e„ - ^ ± a» t a^ ) _ 1 e+ + k z N$e 3 ] 

i,m 

= h Yl I dk ± dk * [&±n+e- + k ± il [ l ) e + + k z IL^e 3 ] , (27) 
i J 

where the operators fl± 3 (k±, k z ) are defined as 

rn 

n w = -iE^-i^. ( 29 ) 

rn 

n? = ( 30 ) 



Notice that the z component of the momentum is diagonal in this basis, just as for plane 
waves, but this is not the case for the other components. Nevertheless, Eq. (127(1 shows that 
Bessel beams with k± ^ may carry linear momentum in the plane perpendicular to the 
propagation direction e 3 . This is the case, for instance, for a field described by two mode 
coherent states such as \ct)i,k ± ,m,k z W)i,k ± ,m±i,k z ,, where as usual, 

a${kz,k±)\a) i;K = a\a) itK . (31) 

From Noether theorem and the isotropy of space, the following definition of the field 
angular momentum in a volume V and around a point r is obtained : 

J(r„) = ^ y (r - r ) x [E(r, t) x B(r, t))dV (32) 
= J(0) - r x P. (33) 

Due to the presence of the position vector r, this integral may diverge if taken over the 
whole space; this may be the case even if the fields are localized and the integrated energy 
and linear momentum are finite. 

Using Maxwell equations, the total angular momentum can also be written as jjj: 

J(r ) = / E i [(r-r )xV]A i dV+-^- [ E x A dV 
Aire J v Aire J v 

- JL jf E[(r-r„) x A] • ds, (34) 

where summation over repeated indices is implicit and S is the surface enclosing V. The 
first integral involves the differential operator (r — ro) x V, which is usually associated to 
the orbital angular momentum; thus, it is customary to identify 

L(r ) = ! Ei[{v- t ) x V]A, dV, (35) 
Aire J v 



'V 

wit h the OAM of the Hold (51. On the other hand 

1 



E x A dV (36) 



v 



Aire 

is independent of the choice of origin and is identified with the spin of the field. 

However, one should be careful with the above identification of spin and orbital terms 
because they depend on the chosen gauge, whereas physically observable quantities should 
not. This difficulty is commonly avoided using the transverse gauge, V ■ A = 0; if each 
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mode of the EM field has a well defined frequency, that is E w (r, t) = dt Eo(r)e~ tuJt+ ^, then, 
in this gauge, A u (r,t) = (—i/u)~E ul (r,t) for each mode, and L and S can be written in an 
apparently gauge independent form. Once the EM field is quantized, the results obtained 
in the transverse gauge turn out to be consistent with the expected values ±h of the spin 
in plane and spherical symmetries, but this is not the case for other gauge selections j^lsj]. 
The former consistency is also obtained for the angular momentum flux in the transverse 
gauge Q- 

Notice also that, in general, the intrinsic angular momentum of a massless particle cannot 



be defined in an unambiguous way. 



(see, e. g., the discussion in Ref. 



in his classical experiment 
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nstead, the relevant dynamical variable is the helicity 
2 If ) and it is actually this quantity that Beth measured 



Finally, it is important to notice that the integral associated to L in Eq. (j34j) is well 
defined only if the electromagnetic field vanishes faster than r~ 2 . Van Enk and Nienhuis []| 
studied the consequences of quantizing the electromagnetic field in terms of creation and 
annihilation operators related to such localized electromagnetic modes: they have shown 
that even under this boundary condition the corresponding operators L and S cannot be 
identified with angular momentum operators because they satisfy a closed but different 
algebra. However, if the electric and magnetic fields are written using a basis formed by 
non localized modes, there is no natural separation of spin and orbital momentum since the 
integrals are not well defined. 

This kind of difficulties is manifest for Bessel beams. One possible way to overcome the 
problem is to impose boundary conditions on a given surface (a cylinder in this case), but 
such a restriction would break isotropy j^J. The other possibility is to take Eqs. (|3l))) and 
f!36|) as definitions of angular and spin operators, and to carry on the calculations in order 
to study the properties of the resulting operators. We will use the latter approach in the 
following. 

The result for the "orbital" angular momentum quantum operator using the basis of 
Bessel modes turns out to be: 



L(0) = ftj^j dk ± dk z 



■ Kz l \" (*)t - (i) 



~ V (m "^ )a - ta - ie+ + m ^ )e3 



) a m-l a m e - 

(37) 
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i 

where the operators A ± 3 (/c^, are given by 

A+ = iE( m -^)ft( fc i. fc *) a 2( fc i. fc .). ( 38 ) 

m 

A« = -^(m-^^.yaL^,^, (39) 
A« = 5>A>»(fc ± A). (40) 

m 

The above relations have a more complex structure than the one obtained for spherical 
vectors (see Appendix B), but this is to be expected since the latter are explicitly constructed 
to describe the orbital angular momentum. 

Notice also that L z (0) is invariant under Lorentz transformations along the z axis as 
expected; it can be interpreted as an intrinsic operator since it does not depend explicitly 
on k±. On the other hand, L X)J/ (0) are highly dependent on quantum numbers {k±,m, k z }; 
moreover, if we define L± = L x ± iL y , then L + (L-) acts as a lowering (rising) operator 
that changes m — > m — 1 (to — > to + 1). 

For the helicity operator S, we find: 

m 

+ k ± {d^l^ - c^a^K + ik,{a$*ag? - a^a^)e 3 ] (41) 

= h dk ± dk z - [k ± £ + e_ + k ± t_e + + A; 2 S 3 e 3 ] , (42) 
J w 

where the operators S± i3 (/c_|_, k z ) are defined as 

S + = ^E^^-i-^^-i)' ( 43 ) 

m 

S- = ^E^-i^-^^), (44) 

Ss = <E( a S )ta 2 ) - a S )a S )t )- ( 45 ) 
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IV. ALGEBRAIC PROPERTIES OF THE DYNAMICAL OPERATORS. 



The basic dynamical quantities can in general be identified by their algebraic properties. 
Thus, for instance, the components of the linear momentum operator must commute among 
themselves; a direct calculation shows that this is indeed the case: the operators U n defined 
in Eq. (J27)l do commute, [Hi , Hj] = 0, and therefore [Pi, Pj] = as expected. 

However, the components of L and S do not satisfy the commutation relations of angular 
momentum. In fact, it can be seen that: 

[L + ,L 3 ] = hL + , (46) 
[L_,L 3 ] = -hL_, (47) 
[L+,L-] = 2h2 YlJ dk ± dk z ^A 3 . (48) 

On the other hand, all the components of the operator S commute among themselves: 
[Ej , Sj] = 0, so that 

[Si, Sj] = 0, (49) 
and it can also be shown that they commute with the momentum operator: 

[P,S,} = 0. (50) 

These properties are compatible with the identification of S as an helicity operator. 

Furthermore, L 3 commutes with the z-component of the linear momentum operator P, 

[L 3 ,P 3 ]=0. (51) 

while 

[L 3 ,P_] = hP_, (52) 

[L + ,P.] = HP 3 , (53) 

[L + ,P 3 ] = 0, (54) 

[L + ,P + ] = J dk±dk z k z a ( £l 1 a m+1 , (55) 

i,m 

and therefore the algebra of these operators does not close. 
Finally, S3 commutes with L, 

[S 3 ,L] = 0, (56) 
11 



while 

[S+,L+] = -HS Z (57) 

/rk 
dk ± dk z —±Z + (58) 

[S+,LJ\ = -ih 2 J dk ± dk z < ^j-{a^_ 1 a^l 1 - a^a^). (59) 

Summing up, the components of the momentum operator P commute among themselves, 
as it should be, but the algebra they generate with the other two operators L(0) and S is 
not the standard one for the translation and rotation group. This is not unexpected since, 
according to our previous discussion, there is an ambiguity with the decomposition of the 
total angular momentum into spin and orbital parts; and moreover, the "spin" is rather the 
helicity. 

The polarization state of a plane wave with propagation vector k = k 3 e z is completely 
characterized by its Stokes parameters. Their quantum counterparts jl^ are given by the 
operators 

a 1 = d {x)] d {v) +d {v)] d {x \ (60) 
a 2 = i(d^d^ -a^d^), (61) 
a 3 = d^d {x) -a^aW, (62) 
a = flWtflW + aMtaW, (63) 

where the indices x and y refer to linearly polarized plane waves in the corresponding direc- 
tions. The operators {ax, cr 2 , 03} satisfy the algebra of the rotation group: [0*, <jj\ = 2i€ijk<Jk 
up to a factor 2. One can readily extend these definitions to Bessel beams identifying the 
indices with the TE and TM superscripts. Clearly, 02 is the elementary operator appearing 
in S3. Thus, measurements of a for Bessel beams should yield important information about 
their polarization states, just as in the case of plane waves. 

Since £ , P3, L 3 (0),and S3 commute among themselves, they can be simultaneously diag- 
onalized. This can be done by introducing the operators, 

aS>=:-^(ag>±iag>), (64) 

which corresponds to a new basis 



A (±) - J_ ( A ( TA/ ) + i A ( Ti? ) 



■ I m 
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for the fundamental modes. 

At this stage, it is important to compare our results with an alternative selection of basis 
modes that appears in the literature [a 0| • Namely, the following modes 

i (h i 



A%\r,t;k x ,k z ) = A 



2\k z 
i (k\ 



A^{x,t\k u k z ) = A[ c) e+i) m - ^ (jr)VWi e 3 



(65) 
(66) 



where ip m (r, t; k±, k z ) = J m (fcj_p) exp{— iut + ik z z + imtfi}. They are considered to be the 
analogues of right (TV) and left (£) polarized plane wave modes |g,|23j. Their superpositions 
Am ± Am define linearly polarized modes, and they can be written as linear combinations 
of elementary TE and TM modes: 



A {n) 



-4, 



m—l 



OJ I 



A(C) _ AC)f( A (TM)_ {TE) \ 
V OJ ' 



(67) 
(68) 



Within the quantization scheme, this change of basis corresponds to the following definition 
of the annihilation operators: 

1 /./ii .ck 7 



t m+l 



,(£) 

'm-l 



y/1 + (Ck z /u)< 
1 



(2) 



a$ + % — -a. 

OJ 



OJ 



Now, the point is that, although the helicity operator S3 is diagonal in this basis: 

1 + (oj/ck z ) 2 



(69) 
(70) 



S 3 = y dk±dk z 



(71) 



the operators £, P3 and L3 are not diagonal. This can be seen from the fact that: 



+ [1 - (c/cA;,) 2 ] (aSU^x - }> ( 72 ) 

as follows with some straightforward algebra. It should be noticed that it is only in the 
paraxial approximation, k z ~ oj/c, that the second term in this last equation, which is non 
diagonal, does vanish. 
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V. DISCUSSION AND CONCLUSIONS 



Let us summarize the main results obtained with the quantization of Bessel beams. The 
proper values of the set of observables {£, P 3 , L 3 (0), S3} define the possible quantum numbers 
that characterize the Bessel photons: {id,hk z ,mh,±hk z c/id}; their physical interpretations 
are clear: for instance, S3 is the helicity operator. We have also analyzed the role of all the 
dynamical operators appearing within the quantization scheme. It turned out that the three 
components of the orbital angular momentum {L + , L_, L 3 } do not satisfy a closed algebra, 
despite the fact that L 3 is related to a spatial rotation around the z axis; in fact, this algebra 
is not the same as obtained for localized fields in Ref. |7[. This is the price we have to pay 
for not taking the surface terms in Eq. (J34|) into consideration. 

Now, the algebra of the local operators E and B is independent of the gauge and the 
basis set. Accordingly, global bilinear operators of the electromagnetic fields such as P and 
S have commutation relations among all their components that are also independent of the 
basis set; this is guaranteed by the normalization condition that each photon carries an 
energy hid. However, the global operators J and L are defined in terms not only of E and B, 
but also of the position vector r; this term induces a strong dependence of J and L on the 
boundary conditions satisfied by E and B. This fact is illustrated in Appendix B, where the 
equivalent L(0) operator is given in terms of the spherical vector basis, and it is shown that 
it does satisfy the standard algebra. In any case, the algebraic properties of the dynamical 
operators and their commutation relations have physical consequences because they imply, 
for instance, specific uncertainty relations that could be verified experimentally. 

It is also worth mentioning that all the dynamical operators we have studied in this paper 
correspond to global observable quantities. A further analysis of local dynamical quantities, 
such as the tensor describing the angular momentum flux, could elucidate the difference 
between "spin" and "orbital" angular momentum. In fact, it was shown by Barnett [2^ 
that in the classical case, there is a natural separation into spin and orbital parts for the z 
component of this flux, M zz . However, a quantum description should also include the full 
commutation relations of the appropriate separated parts of this tensor. This is particularly 
relevant in the light of recent experiments measuring the rates of spin and orbital rotation 
of trapped particles at different distances from the beam axis jl4|. 

It is now well established that Bessel beams induce the rotation of microparticles trapped 
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in an optical tweezers 
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, 3|. The experiments described by O'Neil et al. [15j use Laguerre- 
Gaussian waves that are circularly polarized in the sense of our Eq. (|66p. When the beams 
are converted into linearly polarized waves by a birefringent trapped particle, the parti- 
cle spins around its own axis in a direction determined by the handedness of the circular 
polarization, while small particles trapped off the beam axis rotate around that axis in a 



direction determined by the handedness of the helical phase fronts |l5j. Now, according to 
our results, a similar experiment with Bessel beams that are superpositions of elementary 
TE and TM modes, A(K) = A^™\K) ± iA^ TE \K), should also induce the spinning of a 
trapped particle around its axis; but, since there is a relation S3 = ±k z c/u for each photon, 
the angular momentum should exhibit a linear dependence on k z for a fixed beam intensity. 
This prediction could be tested experimentally. 

In a future publication, we will investigate the quantum electrodynamic interaction of 
atoms with Bessel beams using the formalism developed in this paper. Particular emphasis 
will be given to further clarifying the role of spin and orbital angular momentum of light. A 
detailed analysis of this interaction should explain why the spontaneous emission of Bessel 
photons by atoms is a strongly inhibited process, as experiments have shown so far. 
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Appendix A. Some useful equations. 

The following formulas are used in order to perform the integrations of terms involving 
Bessel functions. They can be easily obtained from the Hankel transform and anti-transform 
formulas (see, e. g., |2J|)- Namely : 

f°° 1 1 

J J m (kp)J m (k'p)pdp= -5(k- k'), (73) 



and 



f°° 1 1 

J J m (kp)J' m {k'p)p 2 dp = --5\k-k')--5{k-k'). (74) 
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From these last expressions and using the standard recurrence relations for Bessel functions, 
it also follows that 

"oo ^ 2 



and 



Jo i^ J -(M^(A: / p)+pj;(M^(^p)Jrfp = ^(A:-A;'), (75) 

/°° 1 m 4- 1 

J m (kp)J m+1 (k'p)p 2 dp = -S'(k - k') + ~^-S(k - k'). (76) 

Using the above formulas, we can obtain several typical integrals that are used in Section 
III. Using the shorthand notation M = M m (x M ; k±, k z ), M' = M m /(x M ; k' ± , k' z ), etc., it can 
be shown that for the scalar products: 

J M ■ M'*dV = J N • N'W = (2n) 2 -^^5 m , m ,5(k ± - k' ± )5(k z - k' z ), (77) 

and 

J M • N'W = J N • M'W = 0. (78) 

Also: 

J M ■ M'dV = - J N • N'dV = -(27r) 2 -^5 m ,_ m ,5(A; ± - k' ± )S(k z + k' z )e~ 2 ^\ (79) 



and 



J M • NW = J N • MW = 0. (80) 
Similarly for the vector products: 



- / M x N'W = / N x M'W 



and 



And also 



= (27r) 2 — |-[5 m+ i jm /e_ - 5 m _i jm/ e + ] + 

-^<We 3 }<5(A; ± - k' ± )5(k z - k' z ), (81) 

J M x M'W = y N x N'W = 0. (82) 

y (M x N' — N x M')dV = 0. (83) 
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Define now the operator C± = C x ± iC y , with C = —ir x V. Then: 

d 



C± = e 



±id 



d i d 
dp p d<£ 



P 



dz 



It then follows that: 



M' • (£+M*)dV 



(84) 



i(2vr) 



2 UUJ „i(u-u')t x 

Orn m! - 



k I k?k*.. 



-'m,m'+l 



and 



k±7^ ™ h m— 

ak 7 ok] k i 



M'* • (£ + M)dV 



<S(*l- (85) 



.(2 7 r) 2 -^-e^— ')*(y m ,m'-i 



A? i kyk 



d d k. 

' dk 7 2 dk i fc i 



Also: 



y M' ■ (£ + M)dV = 0, 



S(k± - k' ± )5(k z - k' z ). (86) 



(87) 



using the fact that in this formula all Dirac deltas appear multiplied by their arguments, 
and x8(x) = 0. 

Appendix B. Comparison with plane and spherical vector EM modes. 



Using the formula 



2tt 



irrupt ik^p[cos(f>cos ip^+sin (p sin ip^] 



(88) 



it can be seen that the vectors M and N, given by Eqs. and (|12|). and that determine 
the TE and TM Bessel vector potentials, can also be written in the form: 



M(r, t; re) = (-i) m / dW k - r 5(k z - k' z )5(k± - k' ± )e imipk ——0 k , (89) 



N(r,t; re) = -(-i) m I d^k'^'^diK - k')5(k ± - k')e imipk -^6 k 

uik i 



(90) 



where 



$k = cos 6^ cos y?kGi + cos 6k sin (pke 2 — sin dke 3 



(91) 
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and 

(fk = - sin y?kei + cos y9 k e 2 (92) 

are the spherical unitary vectors associated to the angular coordinate 8^ and tp-^ in the space 
of the propagator vectors k. These expressions show explicitly the transverse nature of the 
electromagnetic Bessel modes. They can be considered expansions of Bessel modes in terms 
of plane waves that, as it is well known, diagonalize the momentum operator. Eqs. (j89H90j) 
permit to evaluate the expressions for Bessel modes in terms of the spherical vectors. 

Spherical vectors form a complete basis for transverse electromagnetic fields in free space. 
They are defined by (see, e. g., 2l| ) 



1 

where 



A SU r ) = 7^ / ^ASjjk)^, (93) 



A2- m (k) = — ^(|k| - w )yW(fi), n = p. (94) 
In these equations, the superscript specifies the electric (E) and magnetic (B) modes, and 

O^") = -7^-r^nY ]m (0n^n), (95) 
J {J + l ) 

YgPiOiM) = fix^'fe^), (96) 

with 

„ - d 1 d 

V A = ^ k — + ^—5 ; 97 

at^k sm t/k 

Fj m (^, yjfl) are the spherical harmonics. When the electromagnetic field is properly quan- 
tized in terms of spherical vectors (SV) the corresponding angular momentum operator 
l/ 5 ^ takes the form 

1(^(0) = ft £ / [±>/5 - m){j + m + I)a^ im+1 a«. w e_ 

i,j,m 

+ m)(j-m + l)a^ m _ 1 a^ jm e + + mN^ m e 3 , (98) 
with the associated number operator: 

A direct calculation shows that, in these case, the standard commutation relations are 
obtained: [L[ SV \ L^ SV '] — it\£ijkL^ V \ 
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iFrom a straightforward calculation it follows that: 

oo rrij 



N fexfczm (r) = J2J2 ^ u ^ k z,rn\u,j,mj)(A^ m .{T) + Ag^.(r)), 

7=1 -m; 



OO «8j 



with 



r l l 2 h 

-47T 2 (-l) m+(m+H)/2 (2) m+ ^(|k| - w)<J 



'■m.m 



^ k z u 1 / 2 



and 

v{k u k z ,m\uj,j,m j ) 



(2j + l)(j-\m 


)! 9 


47r(j + m )! 


<9/c z 


^ rf 3 rM fc±fezm • A 


(E) 




A (M) 

uijrrij 



= 47T 2 (-l) m+ ( m+ H)/2(,)™+^(| k | _^)^ 



X 



(2j + l)(j - H) ! imu \ m \,ck g 



(100) 



M k±kzm (r) = J2H /^(^,^,m;^,i,m J )(A^(r)-A^.(r)), (101) 

7 = 1 -raj 



P' m p) (102) 



-Pj m| (^), (103) 



47r(j + |m|)! c/cj_ J u 

where Pj'" 1 are the standard Laguerre polynomials. Thus, as expected, Bessel modes are an 
infinite superposition of spherical waves with different orbital angular momentum j. 
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